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ABSTRACT 


This  paper  presents  an  analytical  study  of  the 
interaction  between  incident  wave  fields  and  a  single 
rib  on  a  fluid-loaded  panel.  The  panel  is  modeled  as 
an  infinite  membrane  (with  frequency  dependent  tension 
to  partially  simulate  the  dispersion  characteristics 
of  a  thin  elastic  plate),  and  the  incident  waves  are 
taken  as  plane  structural  or  acoustic  waves  at  normal 
and  oblique  incidence  on  the  rib.  Our  principal  con¬ 
cern  is  with  the  structural  wave  field  transmitted 
across  the  rib  (in  the  case  of  infinite  mechanical 
Impedance  and  finite  nonlocal  rib  impedance)  though 
the  nonspecular  acoustic  field  scattered  by  the  panel- 
plus-rib  is  also  examined.  Matched  asymptotic  expan¬ 
sions  are  used  to  construct  analytical  approximations 
covering  the  entire  frequency  range,  utilizing  the 
smallness  of  a  "fluid-loading-at-coincidence"  parame¬ 
ter.  The  analytical  results  recover  numerical  results 
obtained  elsewhere  in  appropriate  regions  of  parame¬ 
ter  space,  complement  them  by  providing  simple  approx¬ 
imations  in  those  regions  (typically  of  rapid  variation 
or  of  heavy  fluid  loading)  where  numerical  methods  are 
difficult  to  implement,  and  reveal  the  physical  aspects 
of  fluid  loading  in  effecting  energy  transfer  across 
ribbed  structures. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  DTNSRDC  Contract  Nos.  N68171-78-M-9493 
and  N68171-79-M-8380,  and  DTNSRDC  Job  Order  1-1902-005-01.  Dr.  Crighton 
is  Head  of  the  Department  of  Applied  Mathematical  Studies,  University  of 
Leeds,  Leeds,  England. 

1.  INTRODUCTION 

This  report  attempts  to  give  a  comprehensive  analytical  description 
of  acoustic  and  vibration  phenomena  associated  with  the  interaction 
between  plane  acoustic  or  structural  waves  and  a  single  mechanical  line 
constraint  on  a  fluid-loaded  surface.  Specifically,  the  phenomena 
analyzed  concern  the  transmission  of  free  subsonic  surface  wave  energy 


across  the  constraint  (the  "rib"),  the  generation  of  an  acoustic  field  in 
the  transmission  process  (i.e.,  the  scattering  of  vibration  energy  into 
sound),  and  the  non-specular  field  generated  when  a  plane  propagating 
acoustic  wave  is  incident  upon,  and  reflected  by,  the  surface.  Studies 
of  these  processes  have,  of  course,  been  made  before;  for  the  transmission 
coefficient  see  [1,2],  while  the  scattered  acoustic  field  has  recently 
been  examined  in  [ 3 ]  and  [ 4 ] . 

None  of  these  studies  has,  however,  been  able  to  cover  the  whole 
frequency  range  of  interest  because  they  have  not  exploited  the  fact  that 
an  appropriately  defined  fluid  loading  parameter  £  is  invariably  small. 

This  smallness  can  cause  difficulties  in  numerical  calculations,  as 
happened  for  example,  in  [2]  where  it  was  not  possible  to  evaluate  the 
free  wave  energy  transmission  at  low  frequencies  and  at  frequencies 
"close  to  coincidence."  An  explicit  recognition  of  the  smallness  of  £ 
can,  nonetheless,  be  turned  to  advantage  with  the  aid  of  modern  perturba¬ 
tion  theory.  This  can  be  used  to  construct  a  sequence  of  analytical 
approximations  describing  rib-associated  phenomena,  and  to  ensure  that  the 
approximations  blend  smoothly  into  one  another  so  as  to  cover  the  entire 
range  of  frequencies. 

A  first  step  in  this  direction  was  taken  in  [5]  where  approximations 
for  small  £,  covering  the  entire  frequency  range,  were  derived  for  the 
line  transfer  admittance  A^  and  the  line  drive  admittance  Aq  of  a  fluid- 
loaded  surface.  The  surface  considered  was  a  membrane  with  frequency- 
dependent  tension  to  simulate  the  dispersive  wave  characteristics  of  a 
thin  elastic  plate.  Naturally,  of  course,  such  a  device  fails  to  reproduce 
the  near-field  behavior  associated  with  shear  rather  than  tension,  but 
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it  was  still  felt  that  the  membrane  model  could  give  useful  results — in 
particular,  in  the  form  of  physical  Insight — for  the  role  played  by  fluid 
loading  in  the  interaction  between  structural  waves  and  surface  disconti¬ 
nuities.  This  role  was  discussed  in  [2];  the  essential  point  is  that 
fluid  loading  is  able  to  "mend"  a  surface  discontinuity  in  a  way  that  in 
some  frequency  regimes  is  more  or  less  complete,  so  that  a  structural  wave 
can  experience  total  transmission  across  the  rib — even  if  the  rib  has 
infinite  mechanical  impedance  (so  that  total  reflection  would  occur  in  the 
absence  of  fluid  loading) . 

In  the  light  of  these  remarks  it  is  then  appropriate  to  consider 
whether  the  transmission  might  not  be  reduced  in  the  presence  of  fluid 
loading  if  a  suitable  finite  mechanical  impedance  of  the  rib  were 
utilized.  This  will  be  considered  in  detail  in  Section  4  below.  It 
is  also  appropriate  to  inquire  whether  studies  of  structural  waves 
incident  normally  on  a  rib  adequately  cover  applications  in  which  free 
structural  waves  are  incident  on  a  rib  from  all  possible  directions. 

The  problem  of  oblique  incidence  can  be  readily  related  to  that  for  normal 
incidence,  as  explained  in  Section  5  below,  where  the  requisite  trans¬ 
formations  are  given.  It  turns  out,  however,  that  when  the  transforma¬ 
tions  are  applied  to  the  asymptotic  expressions,  the  transformed  expres¬ 
sions  are  not  uniformly  valid  for  certain  particular  directions  of 
Incidence.  Separate  studies  have  to  be  made  for  those  directions,  and 
the  results  are  briefly  reported  in  Section  5 — briefly,  because  it  turns 
out  that  no  new  phenomena  arise  in  the  oblique  incidence  case.  In 
particular,  in  the  frequency  range  for  which  the  "mending  mechanism"  of 
fluid  loading  is  inefficient  for  normally  incident  waves,  it  remains 
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comparably  Inefficient  for  all  obliquely  incident  waves.  Thus,  there 
appears  to  be  no  significant  filtering  (i.e.,  selective  transmission)  of 
waves  incident  from  preferred  directions,  at  least  when  the  constraining 
rib  has  infinite  mechanical  impedance. 

One  is  then  again  prompted  to  consider  the  effect  of  finite  rib 
impedance,  and  in  the  case  of  oblique  incidence  a  local  reaction  impedance 
is  inappropriate.  The  rib  is  therefore  represented  as  a  beam  (or  rather, 
as  the  plate  has  been  replaced  by  a  membrane,  by  a  string)  with  its  own 
nonlocal  reaction.  At  each  frequency  a),  a  certain  direction  of  incidence 
corresponds  to  resonance  of  the  rib,  the  free  wavenumber  of  the  rib  being 
matched  by  the  incident  wavenumber  parallel  to  the  rib,  and  it  is  shown 
that  there  is  perfect  transmission  of  such  oblique  waves  regardless  of 
fluid  loading  effects.  In  particular,  there  appears  to  be  no  way  of  using 
fluid  loading  effects  to  significantly  inhibit  this  selective  transmission 
of  oblique  waves  which  cause  the  rib  to  resonate. 

In  all  this,  the  membrane  model  of  the  surface  will  be  retained. 
Extension  to  the  case  of  a  thin  elastic  plate  may  be  a  useful  development 
in  due  course,  but  a  more  urgent  need  at  present  is  to  extend  these 
studies  of  plane  wave  incidence  on  a  single  rib  to  cases  where  the  incident 
energy  is  generated  by,  say,  a  force  at  a  finite  (small  or  large)  distance 
from  a  rib,  and  where  there  are  several  regularly  or  irregularly  spaced 
ribs.  These  extensions  will  form  Part  II  of  this  report.  A  completely 
general  formulation  covering  all  these  aspects  (arbitrary  numbers  and 
locations  of  ribs,  arbitrary  driving  excitation,  arbitrary  (homogeneous) 
surface  response  operators,  etc.)  is  being  given  under  separate  cover  [6J; 
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the  point  of  Parts  I  and  2  of  the  present  work  is  to  give  a  quantitative 
assessment,  in  simple  model  configurations,  of  the  various  coupling  mech¬ 
anisms  implicit  in  the  general  formulation. 

In  the  following  sections,  a  surface  or  panel  of  the  membrane  type 
occupies  the  plane  y  =  0,  with  fluid  of  density  Pq,  and  sound  speed  co>  in 
y  >  0,  and  a  vacuum  in  y  <  0  (with  trivial  modifications  to  allow  for  the 
presence  of  identical  fluid  in  y  <  0).  A  single  "rib,"  whose  action  on 
the  panel  may  be  represented  at  each  frequency  ii)  by  a  scalar  admittance, 
lies  along  the  line  y  =  0,  x  =  0,  and  cases  of  normal  incidence  involve 
waves  propagating  in  the  x-direction  only.  For  oblique  incidence,  the 
rib  admittance  is  allowed  to  depend  upon  the  wavenumber  £  in  the 
z-direction,  as  well  as  on  the  frequency.  The  time  factor  is  taken  as 
exp(-i^t),  w  >  0,  and  is  omitted.  Extensive  use  will  be  made  of  methods 
and  specific  results  reported  in  [5];  the  readers  familiarity  with  [5] 
will  be  assumed. 

2.  FREE  WAVE  TRANSMISSION:  NORMAL  INCIDENCE 

A  subsonic  free  surface  wave  with  velocity  expdk^x)  is  incident 
from  x  =  -  00  on  a  single  rib  at  x  =  0.  Write  the  total  velocity  of  the 
panel  as 


vtotal(x)  =  Vo  exp(ik!x)  +  vr<x> 


and  define  the  rib  admittance  Af  so  that  the  force  into  the  fluid  with 
which  the  rib  acts  on  the  panel  is  given  by 


v  ,  (x 
total 


0) 


A  F  . 
r 


5 


The  single  scalar  completely  defines  the  effect  of  a  line  mechanical 

constraint  on  a  membrane  type  of  surface.  Next,  relate  v  , (0)  to  the 

total 

amplitude  of  the  incident  wave  through  the  introduction  of  the  line 
drive  admittance  Aq,  giving  (by  definition) 

v  (0)  =  A  F 
r  o 


so  that 


F  = 


A  -A 
r  o 


This  rib-induced  force  then  itself  produces  a  subsonic  surface  wave  field 
transmitted  to  x  =  ±  and  given  by 


vr(|x|  -+■  °°)  ~  A^F  exp ( ik  | x | ) 


Aot  being  the  line  transfer  admittance  to  large  distances.  The  total 
surface  wave  field  at  infinity  to  the  right  of  the  rib  is  thus 


vtotal 


(x 


+  co) 


A  -A  , 
r  o 


V 

o 


expCik^x) 


so  that  a  transmission  coefficient  may  be  defined  as 


T  =  1  + 


A  -A 
r  o 


(2.1) 


which  is  precisely  Eq.  (16)  of  [2]  in  a  different  notation. 


Now  asymptotic  approximations  have  been  given  in  [5]  for  the  admit¬ 
tances  and  AQ,  utilizing  the  fact  that  an  appropriately  defined  fluid¬ 
loading  parameter  is  small  in  most  applications.  Clearly,  these  approxi¬ 
mations  can  be  used  to  evaluate  T  for  any  prescribed  value  of  the  complex 
rib  admittance  A^.  Here,  in  order  to  emphasize  the  role  played  by  fluid¬ 
loading,  we  evaluate  the  transmission  for  a  rib  of  infinite  mechanical 


impedance 


of  the  fluid  coupling.  We  use  the  dimensionless  parameters  of  [51,  namely 


ft  =  oj/oi  ,  e  *  p  c  /rcuti  , 

g  O  o  g 


the  specific  mass  and  coincidence  frequency,  m  and  10^,  respectively, 
being  assignable  quantities.  The  transmission  coefficient  is  then 


T  =  1 


2  e  B°° 


(2.2) 


where  the  dimensionless  admittances  B  ,  B  were  evaluated  in  [5],  The 

oo'  o 

transmission  coefficient  will  now  be  evaluated  for  various  regimes  in  the 
dimensionless  frequency  domain  ft . 

If  e  «  1  and  0  <  ft  <  1  (but  ft  not  close  to  0  or  1)  we  find 


(2.3) 


where  D  ■  1  -  ft.  The  0(e2)  term  can  calso  be  calculated  from  results 
given  in  [5],  A  numerical  evaluation  of  |T^/e|  is  given  in  Table  1,  and 


TABLE  I  -  VALUES  OF  |Tfi/e|  CALCULATED  FROM  EQUATION  (2.3) 


ft 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

lyel 

4.59 

3.28 

2.91 

2.88 

3.10 

3.63 

4.72 

7.52 

18.34 

the  results  agree  with  the  numerical  plot  given  in  [2,  Figure  2]  by 
Maidanik,  et  al,  except  when  ft  is  small  or  close  to  unity.  The  existence 
of  a  minimum  value  of  |T^/e|  close  to  0  *  0.4  should  be  no tea.  Note  also 
that  Eq.  (2.3)  predicts  infinite  transmission  as  ft  0,  with 


T 


ft 


21 

IT 


and  infinite  transmission  4s  ft  -*■  1  -  ,  with 


'n 


e 

2(1  -n)3/tL 


These  predictions  are  unacceptable  and  clearly  indicate  that  the  approxi¬ 
mations  leading  to  Eq.  (2.3)  are  not  uniformly  valid.  We  find,  in  fact, 
that  complementary  "matching"  expansions  are  needed  when  1 1— Q j  =  0(e  n) 
and  when  ft  =  0(e). 

Note  finally  that  Eq.  (2.3)  implies  that  only  a  very  small  fraction, 
O'*  of  the  incident  energy  is  transmitted  when  0  <  ft  <  1. 

We  turn  next  to  ft  >  1  (but  ft  not  close  to  1)  and  find  that 


ft 


1 


j-l si _ \ 

\  ft3/iz(ft-iW 


+  0(e3). 


(2.4) 


Again,  the  0(e2)  term  can  be  obtained  from  [5],  but  serves  only  to  show 

that  Eq.  (2.4)  remains  valid  for  arbitrarily  large  values  of  (2.  Thus  the 

transmission  is  virtually  total  at  all  frequencies  above  coincidence,  a 

surprising  result  at  first  sight,  and  one  not  in  agreement  with  the 

numerical  evaluation  given  in  Figure  2  of  Maidanik  et  al  [2]. 

* 

Substantial  transmission  was  found  there  only  in  the  immediate  vicinity 

of  £2  =  1,  the  transmission  decreasing  rapidly  as  £2  increases  beyond  1. 

Recall,  however,  that  the  results  in  Eqs.  (2.3)  and  (2.4)  apply  only  to 

the  transmission  of  the  subsonic  surface  wave  whether  £2  <  1  or  £2  >  1. 

When  £2  >  1,  the  wavenumber  is  only  slightly  greater  (by  0(t  ))  than  the 

acoustic  wavenumber  k  ,  and  therefore  the  acoustic  impedance  p  w/(k  -k  )^ 

o  o  1  o 

is  very  high.  This  is  reflected  in  a  very  low  value  of  when  £2  >  1, 
whereas  the  drive  admittance  does  not  change  so  dramatically  as  £2  increases 
through  unity;  and  the  effect  of  the  rib,  represented  in  T  by  A  /A  ,  is 

GO  O 

therefore  small  when  £2  >  1.  In  other  words,  because  the  wavenumber  k  is 
only  marginally  subsonic,  the  rib  can  have  little  effect  on  it,  and  fluid 
loading  completely  "mends"  (cf.  [ 2  J )  the  surface  discontinuity  presented 
by  the  rib. 

Now  in  the  calculation  of  Maidanik  et  al  [2],  a  Fourier  integral 
for  the  transmission  at  a  point  x^  >  0  of  the  total  field  generated  by  a 
prescribed  force  at  <  0  was  evaluated  numerically  for  large,  but  finite, 
values  of  xq  and  x  ,  no  attempt  being  made  to  separate  the  field  out 

*Private  communication  from  G.  Maidanik;  in  the  published  paper  the 
curve  for  £2  >  1  was  inadvertently  omitted. 


into  its  various  components  and  to  examine  the  separate  transmission  of 

each  of  these.  When  ft  <  1,  the  dominant  mode  generated  and  transmitted  is, 

of  course,  the  subsonic  free  mode  so  that  our  results  and  those  of  [2] 

would  be  expected  to  agree.  When  ft  <  1,  however,  the  dominant  wave  in 

the  surface,  at  any  rate  for  moderate  distances  x  ,  x, ,  is  the  "leaky  wave 
”  o  1 

field"  (see  [7]),  rather  than  the  subsonic  wave.  Ultimately,  all  the  leaky 
wave  energy  is  lost  to  the  radiation  field,  and  we  are  left  (as  xo*xj  °°) 
with  just  the  subsonic  wave,  but  the  distances  must  satisfy  k^  x  »  e  1 
for  this  to  be  the  case.  One  is  led  to  believe,  therefore,  that  the  calcu¬ 
lations  of  Maidanik  et  al  actually  refer  to  the  transmission  of  the  leaky 
wave  across  the  rib.  To  substantiate  this,  use  the  expression  (4.3)  of 

[5]  for  the  transfer  admittance  to  large  (k  x  »  1)  but  not  too  large 

m 

(k  x  e-1)  distances  of  the  leaky  wave,  and  in  place  of  (2.4)  one  finds 
m 


T (leaky) 


=  -'i/2D&]X 

{l  +  i  -  ?  arc  +  °<e!> 


(2.5) 


This  prediction,  for  | T/e  | ,  agrees  very  well  with  the  results  of  Maidanik 
et  al  [2].  The  transmission,  given  in  Table  2,  is  substantial  only  near 
ft  *  1,  dropping  to  e  e  -’1/2/<2\)  -  0.225  e  e  a r  ft  ®  2  and  decreasing 

like 


e-*i/2 


as  ft  -►  00 . 
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TABLE  2  -  VALUES  OF  | ?^leaky^ /e |  CALCULATED 
FROM  EQUATION  (2.5) 


a 

m 

IS 

m 

IS 

09 

—1 

D 

1.8 

IS 

- ! 

2.0 

|Tleaky/£| 

13.96 

4.44 

2.20 

1.32 

0.87 

0.62 

0.46 

0.35 

0.28 

0.225 

We  argue,  however,  that  it  must  not  be  concluded  that  only  the  lealq 
wave  is  significantly  excited  when  >  1.  That  conclusion  might  be  reached 
from  study  of  a  force-excited  panel  where  the  input  spectrum  is  white  in 
the  wavenumber  domain,  and  where  the  lower  impedance  of  the  leaky  wave 
(compared  with  the  subsonic  wave)  ensures  that  the  latter  is  only  weakly 
excited.  If,  on  the  other  hand,  the  input  spectrum  is  not  white,  but  is 
heavily  concentrated  in  the  subsonic  wavenumber  range  k > kQ,  with  hardly 
any  energy  in  k<kQ,  then  it  is  possible  that  the  subsonic  wave  will  be 
favored  over  the  leaky  wave.  This  may  be  the  case  with  hydrodynamic 
excitations  of  the  system  (turbulent  boundary  layer  pressures  for  example) , 
though  such  excitations  will  generally  have  less  spectral  energy  around 
the  marginally  subsonic  wavenumber  kj  than  at  the  characteristic  hydro- 
dynamic  wavenumber.  Consequently  it  is  essential  that  the  various  wave 
modes  generated  should  be  separated  out  the  the  transmission  established 
separately  for  each.  More  will  be  said  about  this  aspect  (and  about  mode 
conversion,  subsonic  to  leaky  and  vice  versa,  by  the  rib)  in  Part  II  of 
this  paper. 
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Next  the  situation  near  coincidence  is  analyzed,  and  an  approximation 
is  found  which  matches  smoothly  into  the  small  value  *  0(e)  of  (2.3)  for 
ft  <  1  and  into  the  unit  value  of  for  ft  >  1.  The  transmission  refers, 
for  ft  <  1  and  ft  >  1,  to  the  subsonic  free  wave.  The  nonuniformity  is 
characterized  by 


|l-ft|  =  0(e2^) 


and  described  by  approximations  in  which  the  magnified  frequency  variable 


A  - 


1-ft 

e2* 


(2.6) 


is  held  fixed  as  e  +  0.  Using  the  results  of  Section  5  of  [5]  one  finds 


(3p^+2A) 


(2.7) 


plus  a  remainder  term  which  is  0(e  ^3)  if  A  >  A  *  (27/4)  ^  but0(e  ^ )  if 

c 

A  <  A  (including  A  <  0).  Here  p  (A)  is  the  (unique)  real  positive  root 
c  1 

of 


F(q  )  =  q3  +  Aq2  -  1  -  0 
•  o  o 

and  is  one  of  the  roots  given  explicitly  by 


(2.8) 
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3  +  f  co8(l  arc  co,(^7  -  «)} 

J +  f  c°8{?  -  i  «<=  C08(^r-  x)j 

I  +  f  c.s[f  -  i  arc  co8(H-  -  l)J 


(2.9) 


Around  the  critical  value  A^  there  is  a  mild  nonuniformity  at  higher 
order,  but  this  does  not  affect  leading  order  terms. 

Observe  that  is  real  and  0(1)  when  A  =  0(1).  It  can  be  shown  (see 
[5]  for  details)  that  as  A  -*+<»,  expression  (2.7)  decreases  to  the  value 
e/[2(l-ft)  ^ ] ,  matching  imperceptibly  the  asymptotic  behavior  of  as  ft 
increases  towards  1  from  below.  Expression  (2.7)  also  matches  smoothly 
to  the  unit  value  of  for  ft  >  1  (Indeed  a  higher-order  matching  to  the 
two  terms  quoted  in  (2.4)  was  carried  out  in  [5]).  To  see  the  general 
features,  note  that 


—  t 
3A  *A 


Pj (3pt+2A) 


from  which  it  follows  (since  3p^+2Ai<0)  that  increases  monotonlcall^ 
from  0  to  1  as  A  decreases  from  +  00  to  - 


M 

* 


In  summary  then,  the  transmission  coefficient  rises  smoothly  and 
monotonically  from  a  value  0(e)  to  essentially  unity,  the  transition  being 
complete  over  the  range  |l-ft|  *  e^3.  Precise  values  are  easily  computed 
from  (2.7);  in  particular  we  have  the  interesting  result  that 

Tx*i+0(el/3)  (2.10) 

at  X  =  0  (where  the  roots  of  (2.8)  are  p^  ■  1,  p^  *  e2^i/3 ,  p  «  e4iri/3)  ^ 
i.e.,  exactly  at  ft  *  1. 

We  turn  next  to  the  nonuniformity  in  (2.3)  as  ft  ■+■  0.  It  turns  out 
that  there  is  first  a  nonuniformity  when  ft  -  0(e),  then  a  later  one  when 
ft  «  0(e2).  Frequencies  ft  *  0(e)  are  designated  as  "intermediate"  and  a 
scaling  ft  *  eA  is  employed  together  with  expansions  as  £  4  0  for  A  fixed. 
The  transmission  coefficient  in  this  case  turns  out  to  be 


This  expression  matches  the  formula  (2.3)  as  A  becomes  large  and  ft  becomes 

small.  The  transmission  is  larger,  0(e^),  than  T^,  but  expression  (2.11) 

is  still  not  good  down  to  arbitrarily  low  frequencies,  as  it  continues 

to  predict  infinite  transmission  as  A  -*■  0.  Because  of  the  term  involving 

ic 

In  e,  (2.11)  is  not  quite  in  the  form  (£*)  times  (a  universal  function  of 
A  alone),  but  it  is  effectively  in  that  form  as  the  In  £  term  varies  very 
slowly. 
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Next  the  "low-frequency"  region  SI  ■  0(e2)  Is  examined,  taking 
-  e2u>  with  expansions  for  fixed  u  as  c  0.  The  descriptions  are  compli¬ 
cated  and  dealt  with  in  detail  in  [5,  Section  7].  Define 


G(q  )  =  q2 (1+q  )  -  ~ 

O  O  O  0) 


(2.12) 


and  denote  the  roots  of  G(q  )  *  0  by  s  ,  s  ,  s  .  In  all  cases  s  is  used 

to  refer  to  the  unique  real  positive  root  (ui  being  always  positive) ,  while 

s  ,  s  are  real  and  negative  if  5  >  i  *  27/4  and  form  a  conjugate  pair 
2  3  c 

-  a  ±  16,  with  a,  6  >0,  when  ai  <  Hi  .  Define  also 

c 


3  In  s 

y  _ ljlI _ 

,  s  (3s  +2) 
n*l  n  n 


(2.13) 


and  0  ■  arc  tan  (6/ a),  0  <  0  <  n/2.  Then  for  u>  >  £5  one  has  the  follow- 

c 

ing  expression  for  the  transmission  coefficient 


“■[s  + 


ITT 

Sj (3Sj+2). 


+  0(ez) 


(2.15) 


while  when  0  <  w  <  ui 


r  %  (3si+i>  i 

L  (38,-1)’  (3Sj+2) J 

“  f  Jk  w  1 

L  (3s  -1)*  (3s,+2)  s, (3s  +2 LI 


(2.16) 


(it  being  difficult  to  improve  upon  the  o(l)  estimate  of  the  error  when 
to  <  u^).  After  extensive  manipulation  it  is  possible  to  rewrite  both 
(2.15)  and  (2.16)  in  terms  solely  of  the  root  Sj  (which  determines  the 
free  subsonic  wavenumber) ,  using  the  expressions 


In 


(Sj+1) 


2Sj (3s 1+2) 


for  to  >  to  ,  and 
c 


In 


S  = 


_si+1J 


iwt* 


2s  (3s  +2) 

l  l 


1+s 


J 


l+3s. 


L1-3>iJ 


a+s^  +  d-3Sl)% 

(l+Sj)^  -  (l-3Sj)^ 

2Sj (3Sj+2) 


0  =  arc  tan 


3s  -1 

1  I 

si+1 


(2.17) 


(2.18) 


(2.19) 


for  co  <  to  ,  noting  that  3s  -  1  >  0  according  as  to  $  co  so  that  all 
radicals  in  these  expressions  are  real  and  positive.  These  expressions 
look  complicated,  but  can  be  easily  evaluated  using  the  Table  of  values 
of  s^  given  in  reference  [5].  They  show  that  in  the  frequency  range 
ft  ■  0(e2),  an  0(1)  fraction  of  the  incident  energy  is  transmitted  across 
the  rib,  and  fluid  loading  effects  have  again  managed  to  "mend"  the 
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rib  discontinuity  even  though  the  fluid  loading  parameter  e  (which 
actually  measures  fluid  mass  divided  by  surface  mass  at  the  coincidence 
frequency)  is  small. 

Such  examination  of  higher  terms  in  the  "low-frequency  expansion" 
as  it  is  possible  to  carry  out  indicates  that  there  are  no  further  leading 
order  nonuniformities  (in  the  transmission  coefficient,  at  any  rate),  and 
that  (2.16)  should  be  good  down  to  zero  frequency.  If  the  limit  w  0  is 
taken,  it  is  found,  on  reduction,  that 


T~  -*•  i  e-1Ti/3  as  a  -  0  (2.20) 

U)  L 

This  result  is  in  fact  true,  for  sufficiently  low  frequencies,  whatever 
the  value  of  the  parameter  e.  It  shows  just  how  substantial  the  trans¬ 
mission  becomes  in  the  "heavy  fluid  loading"  limit. 

There  is  a  nonuniformity  in  these  results  at  higher  order  in  a  narrow 
region  around  the  value 


u  =  5  -  27/4  (2.21) 

c 

where  one  defines 

a  =  ^(l+e2£)  (2.22) 

and  considers  =  0(1)  as  e  -*■  0.  The  nonuniformity  does  not  affect  the 
leading  order  terms,  however,  and  it  is  found  that 
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1  -  j  In  A 

~  ,  7  21tt 

1  -  -  In  4  +  -y 


+  0(e) 


(2.23) 


with  |T^|  -  0.25  +  0(e).  This  higher  order  nonuniformity  manifests 
itself  in  the  numerical  computations  of  [21,  where  in  Figure  2  one  notices 
a  narrow  plateau,  centered  on  ft  =  (0.26) 2  -  0.068  which,  fore  =0.1  (taken 
in  [2])  is  almost  identical  with  our  ,  and  at  which  |  T J  has  a  local 

maximum  of  0.27  in  good  agreement  with  the  prediction  above  of  0.25.  As 
remarked  in  [5],  one  might  have  suspected  the  curious  little  plateau  of 
arising  from  some  numerical  inaccuracies,  but  the  analysis  here  shows 
that  it  is  indeed  genuine. 

The  overall  pattern  of  free  subsonic  wave  energy  transmission  is 
summarized  schematically  in  Figure  1,  in  which  is  also  included,  for 
ft  >  1,  the  transmission  coefficient  for  leaky  waves.  The  various  approxi¬ 
mations  were  shown  in  [5]  to  overlap  in  the  manner  indicated  in  Figure  1. 
This  Figure  largely  speaks  for  itself,  and  the  only  further  comment  which 
needs  to  be  made  is  that  it  is  possible  to  identify  the  physical  mechan¬ 
isms  which  dominate  the  various  approximations.  Inasmuch  as  it  is 
possible  to  express  the  admittances  entirely  in  terms  of  a  wavenumber  kj 
of  a  free  subsonic  surface  wave,  we  need  only  quote  the  mechanisms 
which  dominate  in  determining  that  wavenumber  (in  contrast  to  [5],  where 
a  slightly  more  elaborate  description  was  given).  Our  model  has  retained 
the  following  five  possible  physical  attributes:  (1)  membrane  stiffness, 
(II)  membrane  mass,  (III)  fluid  pressure  forces,  (IV)  fluid  compressi¬ 
bility,  and  (V)  fluid  mass.  For  0  <  ft  <  1  th«  dominant  processes  are  (I) 
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:  1  -  Schematic  Presentation  of  Transmission  Coefficient 
Magnitude  (for  z  «  1)  as  Function  of  Frequency. 


and  (II);  for  ft  >  1  they  are  (III)  and  (IV)  for  the  subsonic  wave,  but  (I) 
and  (II)  for  the  leaky  wave;  for  |l-ft|  =  0(e  ^)  all  five  processes  are 
significant;  for  ft  =  0(e)  only  (I)  and  (II)  are  important;  for  ft  =  0(ez) 
only  (IV),  fluid  compressibility,  can  be  neglected;  and  when  ft  0  the 
mechanisms  (I),  (III),  and  (V)  determine  the  heavy  loading  limit. 

This  completes  our  study  of  normal  incidence  energy  transmission 
across  a  single  rib;  oblique  incidence  will  be  discussed  in  Section  5. 
We  turn  now  to  acoustic  aspects  of  the  action  of  the  rib. 


3.  THE  ACOUSTIC  FIELD  SCATTERED  BY  A  RIB 
In  this  section  we  deal  with  the  acoustic  fields  generated  when  a 
free  subsonic  wave  is  incident  normally  upon  the  rib,  and  when  a  plane 
acoustic  wave  is  incident  upon  the  panel  with  the  single  rib,  the 
latter  being  determined  by  the  so-called  "non-specular  reflection 
coefficient."  Both  acoustic  fields  are  readily  derived  from  the  field 
generated  by  a  single  prescribed  line  force  F  acting  on  the  panel,  that 
field  being  given  by  the  Fourier  integral 


p(x,y) 


Pp^2  F 
2irT 


■+°°exp{ik|x|  -  (kz-k2)^  y) 

(k2-k2) (k2-k2)^  -  pkz 
mo  m 


(3.1) 


The  notation  is  standard  except  for  U,  equal  to  PQ/in  where  m  is  the 

specific  mass  of  the  panel.  Also  one  has  k^  ■  mu»z/T  with  T  the  membrane 

tension,  and  k  «  w/c  .  Evaluation  of  this  integral  by  steepest  descents 
o  o 


or  stationary  phase  gives  an  acoustic  farfield 


p(r,0) 


p  U)  F 

o _ 

27TT 


(2k  tt/r)^  e^^°r  sin8 

_ o _ 

[  (k2  cos20  -  k2)(-ik  sin0)  -  pk2  ] 
o  mo  m 


(3.2) 


with  |x|  =  r  cos0,  y  =  r  sin0,  0  <  6  <  tt/2.  This  estimate  is  valid  for 

all  0  if  r  is  sufficiently  large.  If  e  «  1  and  0  <  fl  <  1  it  is  probably 

enough  to  take  k  r  »  1,  but  when  ft  >  1  the  estimate  holds  away  from  the 
o 

Mach  directions  0  =  arc  cosft  2  when  k  r  »  1,  and  close  to  those  direc- 

o 

-2 

tions  only  when  k^r  is  very  large,  k^r  >>  e  .  The  reader  is  referred 

to  [7]  for  further  details  of  the  Mach  wave  field  associated  with 

leaky  waves  in  the  surface,  especially  for  a  description  holding  when 

—  2 

kQr  is  not  as  large  as  e 

Here  Eq.  (3.2)  will  be  accepted  as  describing  "the  distant  acoustic 
field"  in  an  appropriate,  known  sense.  Then 


lp|2  «  ^7  «(ft,e,0)  (3.3) 

O 

with 

H(ft,e,0)  =  - - -7-  (3.4) 

[  (ft  cos20  -  l)2  sin2@  +  -jj-] 

and  we  now  state,  very  briefly,  leading  order  approximations  to  the 
directivity  function  H. 

(1)  ft  >  1.  Define  a  Mach  angle  0  =  arc  cos  ft  ^  and  write  8  *  0W  +  ex 

M  M 

when  0  is  close  to  0„.  If  1 9—0., I  is  not  as  small  as  0(e) 

M  M  - 


(ftcos2©  -  l)2 
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(3.5) 


except  when  0  is  small,  0  =  0(e)  =  eip  say,  and  for  this  case 


H  ~  - ^ -  (3.6) 

(fi-1) 2  ip2  +  i 

This  vanishes  for  ^  -*  0,  a  manifestation  of  the  well-known  "Lloyd's 
mirror"  effect  for  grazing  incidence  acoustic  waves. 

If  | 0— ©M |  =  0(e),  so  that  x  =  0(1) ,  then 


H 


sin20 


M 


ez[4XV  cos20m  sin40M  +  ±] 


(3.7) 


—  Z 

whose  0(e  )  largeness  demonstrates  the  intense  field  at  very  great 

—  2 

distances  (k  r  »  e  )  within  0(e)  of  the  Mach  directions  as  compared 
o 

with  the  0(1)  fields  (3.5)  and  (3.6)  elsewhere. 


(2)  0  <  Q  <  1.  Here  the  fluid  loading  term  e2/fi  must  be  retained  only 

for  0  =  0(e),  where  it  leads  again  to  the  form  (3.6).  This  low  angle 
nonuniformity  arises  in  all  cases,  and  will  receive  no  further  comment. 
For  0  larger  than  0(e), 


H  ~  - - -  (3.8) 

(l-ftcos20)2 

which  h as  only  a  weak  dependence  upon  0  and  indicates  fairly  uniform 
directivity. 


(3)  | 1-ftl  =  0(e  ^ ) .  No  special  phenomena  arise  in  this  case.  The 

Mach  angle  0W  is  0(e  ^3),  whereas  the  Lloyd's  mirror  effect  operates  only 
M 

for  small  angles  0  =  0(e),  so  that  the  Mach  wave  field  continues  to  be  as 
described  in  (1)  above,  the  field  for  ft  <  1  as  in  (2)  above. 

(4)  ft  °  0(e).  Here,  writing  ft  =  eA  and  letting  e  -+•  0  one  finds 


H  ~  1  (3.9) 

and  the  field  is  isotropic  (except  for  0  =  0(e)). 

(5)  ft  =  0(e2).  Put  ft  =  e2w  and  let  e  -+•  0.  Then  for  all  angles, 
including  0  *  0(e), 


sin2Q 

sin20  +  - 
0) 


(3.10) 


which  indicates  a  significant  effect  of  fluid  loading  on  the  field  at 
all  angles,  that  effect  leading  to  a  dipole  type  of  behavior. 


H  ~  5  sin20 


(3.11) 


when  S  *  1. 

Most  of  the  results  given  here  are,  of  course,  well-known,  though  it 
may  be  helpful  to  have  them  set  out  as  above,  as  this  makes  it  clear 
just  when  the  different  approximations  hold. 

To  relate  these  results  to  the  field  scattered  when  a  subsonic 
surface  wave  is  incident  upon  a  rib,  all  one  has  to  do  is  to  substitute 
the  expression 
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(3.12) 


F 


V 

o 


A  -  A 
r 


given  in  Section  2,  into  (3.3).  Here  Vq  is  the  amplitude  of  the  incident 

wave,  A  the  admittance  of  the  rib,  and  A  the  drive  point  admittance  of 
r  o 

the  fluid-loaded  surface.  The  dependence  of  Aq  upon  ft  and  £  has  been 
thoroughly  explored  in  [5],  while  A^  is  a  specified  function  of  ft  alone. 
The  substitution  of  (3.12)  into  (3.3)  therefore  changes  only  the  level 
of  the  sound  pressure  in  a  fully  known  way,  and  does  not  change  the 
directivity,  which  is  completely  defined  by  H(ft,e,9). 

As  for  the  non-specular  reflection  problem,  suppose  that  a  plane 
acoustic  wave  with  pressure 


p  exp(ik  x  cos0  -  ik  y  sin9  ) 
o  o  o  o  o 


is  incident  at  an  angle  0  to  the  surface.  Then  it  is  easy  to  show 

o 

(a  general  formulation  for  the  case  of  many  ribs  is  given  in  [6])  that 
the  non-specular  field  scattered  by  the  rib  is  again  given  by  (3.1)  with 


where 


F  = 


ik  sin0  (1-R. ) 
_ o _ o  * 

(A  -A  ) 
r  o 


R*  = 


ift^  sin©  (ftcos2  0  -1)  -  £  \ 

n  rt  ' 


ift"2  sin9  (ftcos2  0  -1)  +  £  / 
0  o 


(3.13) 


(3.14) 


is  the  plane  wave  reflection  coefficient.  This  value  for  F  is  of  the 
same  form  as  (3.12),  with  Vq  in  that  expression  replaced  by  the  velocity 
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p  iuj 
o 


(-ik  sin0  ) 
o  o 


which  would  be  generated  at  x  =  0  by  incident  and  reflected  acoustic  waves 
if  the  rib  were  not  present.  The  directivity  function,  of  incidence  angle 
0q  and  observation  angle  0,  arising  from  (3.2)  with  (3.12)  and  (3.14),  is 
called  the  non-specular  reflection  coefficient.  Its  dependence  upon  0 
has  already  been  dealt  with.  The  dependence  of  the  squared  magnitude  of 
this  coefficient  upon  incidence  angle  is  as 


e2  sin20 

_ o _ _ 

e2  +  SI  sin20  (Si  cos20  -l)2 
o  o 


which  is,  as  expected  from  Reciprocity  requirements,  identical  with  the 
0-dependence.  The  behavior  of  the  non-specular  reflection  coefficient  is 
therefore  already  described  by  the  remarks  above  relating  to  H(ft,e,0). 

4.  EFFECTS  OF  FINITE  RIB  IMPEDANCE 
In  the  free  wave  energy  transmission  problem  considered  in  Section  2, 
there  will  be  substantial  direct  "mechanical"  transmission  across  a  rib 
unless  the  rib  has  very  high  mechanical  impedance.  We  have  seen  that  in 
certain  frequency  regimes  there  is  also  substantial  transmission,  by 
virtue  of  fluid  loading  effects,  even  when  the  rib  has  infinite  impedance. 
The  question  to  be  addressed  here  is  the  following:  "In  those  frequency 
regimes,  is  it  possible  to  reduce  the  transmission  by  balancing  fluid¬ 
loading  effects  against  finite  impedance  effects?" 
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In  terms  of  the  impedance  Z^  «  A^1  of  the  rib,  the  normal  incidence 
transmission  coefficient  is 


1  +  Z  (A  -A  ) 
_ r  00  o 

1  -  Z  A 
r  o 


(4.1) 


in  which  we  write 


Z  =  -R  -  iX 
r 


(4.2) 


for  the  resistive  and  reactive  components,  so  that  for  a  simple  damped 
mass-spring  system 


X  =  -  (K-Mw2) 
u' 


(4.3) 


and  K,  R,  M  are  all  positive  constants. 

Near  coincidence,  i.e.,  for  |l-fl|  =  O(e^),  there  is  significant 
transmission,  and  it  is  found  from  the  admittance  approximations  given 
in  [5]  that  here 


1  +l 

m 

1  2  (  3p,+2X  )  (R+1X) 

/ 

?  k  V> 

IV_  4  1 

1+US)i(R+ix) 


The  quantity  Pj/(3pj+2A)  is  the  transmission  coefficient  for  infinite  rib 
impedance,  and  increases  from  very  small  positive  values  to  unity  as  the 

2a, 

frequency  increases  through  a  range  0(e  '  )  around  the  coincidence  frequency. 
It  is  clear  that  there  can  be  no  choice  of  R,X  which  will  lead  to  very 
substantial  reductions  In  the  value  of  |T|  in  the  range  where  pJ/(3p1+2X) 
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is  close  to  unity.  That  is  really  not  surprising,  for  around  the  coinci¬ 
dence  frequency  the  surface  wave  is  essentially  a  plane  acoustic  wave  in 

the  fluid,  with  very  high  impedance  p  o>/ ( k2-k2)^  in  directions  normal  to 

o  o 

the  surface.  The  rib  will  therefore  not  be  able  to  significantly  affect 
the  propagation  of  that  wave,  whatever  the  rib  impedance. 

Consider,  on  the  other  hand,  very  low  frequencies,  where  fluid  load¬ 
ing  effects  are  dominant  and  where  the  inf inite- impedance  transmission  is 
given  by  (2.20).  Here  one  finds 


T  1  -  3^(^)al1*  <*+«> 

(R+iX> 

from  which  it  can  be  seen  that  the  choices 


(4.5) 


R  =  0, 


(4.6) 


minimize  T,  and  in  fact  make  T  »  0.  The  negative  value  of  X  required 

for  (4.6)  implies  that  the  rib  is  driven  below  its  resonance  frequency; 

the  balance  expressed  by  (4.6)  is  then  essentially  between  the  stiffness- 

like  behavior  of  the  imaginary  part  of  the  surface  drive  point  impedance 

and  the  mass-controlled  imaginary  part  of  the  rib  Impedance.  (Attention 

was  drawn  in  [5]  to  the  surprising  fact  that  In  A  <0  for  the  fluid- 

o 

loaded  membrane  throughout  the  whole  frequency  range.  Implying  that  fluid 
loading  acts  as  an  equivalent  stiffness.)  Equation  (4.6)  (or  its  analog 
for  a  more  realistic  type  of  panel)  might  only  be  capable  of  satisfaction 
over  a  very  narrow  frequency  band,  if  at  all;  there  is,  nonetheless,  a 
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CCfeMfc-*' 


clear  possibility  of  reducing  the  low-frequency  free  wave  transmission  by 
using  a  mass-controlled  rib  impedance  to  at  least  partially  cancel  the 
fluid-loading  Induced  drive  point  stiffness. 


5.  FREE  WAVE  ENERGY  TRANSMISSION;  OBLIQUE  INCIDENCE 
The  rib  lies  along  the  z-axis,  and  we  now  assume  that  the  whole  wave 
field  has  the  dependence  exp(iilz),  so  that,  for  example,  Jl  is  the 
z-wavenumber  of  a  free  surface  wave  obliquely  incident  on  the  rib. 

It  is  then  easy  to  see  that  any  problem  with  this  z-dependence  can 
be  reduced  to  a  purely  two-dimensional  one  in  the  (x,y)  plane,  provided 
that  in  the  two-dimensional  problem  we  make  the  replacements 


and 


k  -*•  K 
m  m 


(k2-*2)^ 

m 


k  ->•  K  *  (k2-£2)^  \ 

o  o  o  ( 


(5.1 


p  -*•  v  =  pk2/K2  I 

mm  I 

Since,  in  the  original  integral  defining  the  surface  response  (Eq.  2.10 

of  [5]),  the  points  kQ  and  k^  must  lie  above  the  path  of  integration,  it 

follows  that  we  must  take  K  *  +  i(£,2-k2)\  K  *  +  i(Jl2-k2)^  when  K  and 

m  m  o  o  m 

Kq  are  not  real.  For  normalized  functions  which  depend  only  on  the  two 
parameters  ft  and  £  (equivalent  to  the  ratios  of  the  three  wavenumbers 
in  (5.1))  the  correspondence  (5.1)  may  be  expressed  as 


where 


(5.2) 


l  -  k  L 
m 


Consider  first  the  wavenumber  in  the  x-direction  of  a  free  surface 
wave  at  frequency  w;  for  that  wavenumber  to  be  real,  1  must  not  exceed 
the  total  wavenumber  magnitude  of  a  free  wave.  In  reference  [5]  approxi¬ 
mations  were  given  for  the  total  wavenumber,  so  that  one  has  the  following 
values  of  the  maximum  permissible  z -wavenumber : 


«  >  1: 


0  <  fi  <  1: 


l  =  k 
max  o 


0* 


£  +  0(e“) 


2tt2 («-l)2 


>} 


max 


+  -  +  0(e2)  | 

“  v.  2JP(i-nr  J 


A  »  |  =  0(1): 


l  -  k  / 1  +  +  0(e) 

max  m  1  2A^ 


Si  *  —  ■  0(1) : 

e2 


A  *  k  /S i^  s. 
max  m  * 


(where  S!  is  the  positive  root  of  sJ  +  s2  *  ^5  ■  0) 


to  -*•  0  : 


SL  ~  k  /Si 

max  m 


Si‘*  - 


(5.3a) 

(5.3b) 

(5.3c) 

(5.3d) 

(5. 3e) 
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A  =  ^ -  =  0(1):  £  -  k  (l+%e2/j  p.-HUe**)}  (5.3f) 

_/3  max  m  1 


(where  pj  is  the  positive  root  of  p3  +  Ap2  -1=0). 

One  sees  then  from  (5.3)  that  for  frequencies  0(e2)  there  is  a 

change,  by  an  0(1)  factor,  in  the  maximum  permissible  £,  from  its  value  k 

m 

in  the  absence  of  fluid  loading,  and  that  as  -*■  0  this  factor  tends 

(slowly)  to  infinity.  In  this  "heavy  fluid  loading"  limit,  all  waves  can 

propagate  in  the  x-direction,  no  matter  how  large  their  z-wavenumber.  In 

all  other  frequency  ranges  below  coincidence  the  fluid  loading  correction 

to  the  vacuum  value  k  is  small,  and  corresponds  to  an  increase  associated 

m 

with  the  addition  of  fluid  mass  to  surface  mass.  For  ft  >  1,  a  wave 
propagates  in  the  x-direction  with  no  attenuation  if  (5.3a)  is  satisfied, 
and  the  wave  is  subsonic.  If,  however,  the  wave  is  an  oblique  subsonic 
leaky  wave,  its  total  wavenumber  is 


+ 


le 

2ft^(n-n^ 


+  Q(e2)  » 


and  therefore  it  will  propagate  in  the  x-direction  with  only  slow 

attenuation,  0(exp(-Gx)),  if  £  <  k  ,  but  with  strong  attenuation  if 

m 

£  >  k  .  Thus  we  can  take 
m 


ft  >  1: 


£ 


max 


k 

m 


(5.3g) 


as  the  condition  for  a  leaky  wave  to  "propagate"  in  the  x-direction. 
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Results  for  the  transmission  coefficient  acrosss  a  rib  of  infinite 


impedance  were  presented  in  Section  2,  as  functions  of  e  and  0.  It  is 
straightforward  then  to  use  (5.2)  to  immediately  produce  results  for 
oblique  incidence,  though  the  results  are  algebraically  complicated. 
Suppose  first  that  Q  >  1;  then  if  L2  <  1,  ft*  >  1  and  therefore 
T  *  1  +  0(e2)  from  (2.4),  while  if  L2  >  1,  (2^  <  1  and  therefore  T  *  O(e^) 
from  (2.3).  Since  e*  =  0(e)  unless  L2  is  close  to  £2  or  to  1,  we  can  say 
that  when  1  >  1  an  oblique  wave  with  L2  <  1  is  essentially  totally  trans¬ 
mitted  (as  for  normal  incidence),  but  has  a  transmission  coefficient  of 
only  0(e)  if  L2  >  1. 

Suppose  next  that  0  <  Q  <  1,  with  0  not  close  to  0  or  1.  When  L  is 
close  to  1  a  separate  examination  is  needed  anyway,  so  that  we  can  say, 
from  (5.3b)  that  L2  <  1.  Then  £2*  <  1,  and  Eq.  (2.3)  gives  the  trans¬ 
mission  coefficient  except  when  £2*  is  small,  that  is,  when  L2  is  close 
to  £2.  We  have  made  a  detailed  analysis  of  the  latter  possibility,  using 
a  scaling  L2  ■  £2  -  pe2  and  keeping  p  fixed  as  e  0.  By  direct  calcula¬ 
tion  of  the  case  L2  =  (2  we  confirm  that  expansions  in  the  p-variable  are 
in  fact  valid  down  to  p  *  0.  Throughout  the  region  p  =  0(1)  it  is  found 
that  the  admittances  are  given  by 


A 

o 


U) 


Ik«^?  + 


0(e) 


(5.4a) 


A°°  ■  tIt  +  0(c)  (5-4b) 

■  2D* 

and  hence  that  T  ■  0(e).  Specifically,  when  L2  •  £2  we  find 
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W* 


T'^'(t)[^+1”(t)j+0<v2/K-)  <5'4c) 

with  v/K  =  (e/ft^)(l-L2)_%. 
m 

The  prediction  T  =  0(e)  comes  also  from  use  of  (2.3)  away  from  L2  =  ft, 
so  that  we  conclude  that,  except  perhaps  for  L2  ~  1,  but  includ ing  the  case 
L2  ~  ft,  the  transmission  coefficient  is  small,  of  order  e,  for  both 
normally  and  obliquely  incident  waves  with  frequencies  satisfying  0  <  ft<l. 

As  ft  decreases  further,  the  transmission  coefficient  increases, 

%  2 

becoming  0(e)  when  ft  =  0(e)  and  0(1)  when  ft  =  0(e  ),  just  as  for  normal 

incidence.  The  functional  forms,  of  course,  vary  with  L,  but  this  does 
not  change  the  overall  picture,  and  in  fact  as  ft  -*•  0  the  dependence  on  L 

drops  out,  and  the  transmission  coefficient  is  given  by  (2.20)  for  all 

waves. 

The  case  L  ~  1  remains  to  be  dealt  with.  It  turns  out  that  two 
"inner"  expansions  are  needed  to  cover  the  nonuniformity  here,  reflecting 
the  fact  that  the  basic  approximation  (2.3)  is  known  to  be  nonuniform 
when  e  ~  ft  and  when  e2  ~  ft.  If  e^  ~  ft^  as  L2  -*■  1  one  finds  that 

L2  -  1  =  0(e),  while  if  e2  ~  ftA  as  L2  -*■  1  one  must  have  L2  -  1  =  0(e  ^3) . 

2  2a 

One  therefore  defines  L  =  1  -  ye  '  ,  say,  and  constructs  approximations 
to  the  admittances  with  y  fixed  as  e  4  0.  These  approximations  cannot 
be  gotten  from  approximations  like  (2.3);  indeed  the  point  of  introducing 
the  y-variable  is  to  remedy  the  deficiencies  of  (2.3)  when  Y  is  0(1),  so 
that  the  new  approximations  have  to  be  derived  from  exact  results  given 
in  [5],  and  then  matched  to  the  limiting  behavior  of  (2.3).  The  process 
is  straightforward,  but  lengthy,  and  leads  to  admittances 
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Ao  =  Tk“{  ~i~T  -  ~  -m  +  0(el/3) 
°  Tkm  I  2y^cV3  4y/2^  > 


(5.5) 


Tk2  4f 

in 


Fu^f1  -  ~^  +  °<£2/3>  i 

'2e/3  l  2YW  J 


(5.6) 


from  which  we  find 


T  -  1 


OO  V, 

-  T-  =  °<^ > 


(5.7) 


We  have  not  carried  the  calculation  far  enough  to  see  if,  in  fact,  the 
coefficient  of  the  0(e^3)  term  in  T  is  zero  (which  would  then  lead  to 
T  *  0(e)),  but  at  any  rate  it  is  clear  that  the  transmission  is  still 
small  when  T  =  0(1). 

The  approximations  (5.5)  and  (5.6)  do  not  apply  right  down  to 

2  Vi 

T.  =  1  (i.e.,  Y  =  0),  as  the  ratio  of  consecutive  terms  is  0(e  /y)  and 

not  ur.iformly  small.  To  remedy,  this  one  writes  L2  =  1  -  Kt,  constructs 
approximations  for  tc  =  0(1)  and  matches  them  to  the  asymptotic  behavior 
of  (5.5)  and  (5.6)  as  y  ■+■  0.  The  results  are 


Tk  _  k 

m  2  e  2 


+  0(e^) 


(5.8) 


-4  ^ + ~rV +  o(c"^> . 


(5.9) 


Tk  4e 
m 


T  =  0(e) 


and  these  results  appear  to  apply  right  down  to  <  «  0,  so  that  no  further 


33 


"inner"  expansions  are  needed.  Although  the  dependence  of  T  upon  K  has 
not  been  obtained,  even  to  leading  order,  it  was  not  felt  to  be  worthwhile 
to  pursue  this  matter  further.  The  essential  point  is  that  there  are 

nonuniformities  around  L2  =  £2  and  L2  =  1  (reflected,  for  example,  in  the 

* 

much  larger  values  of  the  admittances  as  L2  -*■  1) ,  but  these  nonuniformities 

affect  the  drive  admittance  A  and  the  transfer  admittance  A  in  the  same 

o  00 

way,  and  still  leave  a  very  small  value  of  the  transmission  coefficient 
(for  a  rib  of  infinite  mechanical  impedance)  both  around  Lz  =  ft  and  around 
L2  =  1. 

If  the  rib  impedance  is  not  infinite  there  are  many  possibilities, 

but  in  that  case  there  will  be  a  substantial  mechanical  transmission  of 

energy,  largely  independent  of  fluid  loading  effects.  What  has  been 

shown  here  is  that,  where  fluid  loading  is  the  only  mechanism  for  energy 

transfer,  no  dramatic  changes  are  associated  with  oblique  incidence;  in 

particular,  waves  with  z-wavenumber  close  to  either  k  or  k  are  not 

o  m 

transmitted  with  significantly  different  efficiency  from  normally  incident 
waves . 

In  the  case  of  finite  rib  impedance,  this  will  normally  depend  on  the 
z-wavenumber  and  on  the  frequency  separately  if  the  rib  is,  for  example, 
of  the  string  or  beam  type  with  nonlocal  reaction.  Thus  in  (4.3)  the 
stiffness  K  could  be  written  as 


Mlo2U) 

o 

where  the  resonance  frequency  is  a  function  of  Z.  For  certain  ranges 
of  driving  frequency  u)9  the  condition  oj  s  w  nay  bo  not  for  a  certain 
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small  range  of  5L,  and  in  that  case,  provided  the  resistive  part  of  the 

rib  admittance  is  then  small  compared  with  A  and  A  ,  it  follows  from  (4.1) 

r  00  o 

that 


T  s  1 


(5.10) 


This  almost  perfect  transmission  is  the  result  of  mechanical  coupling 
across  the  rib,  and  one  may  ask  whether  fluid  loading  effects  can  sig¬ 
nificantly  change  either  the  magnitude  of  the  transmission,  or  the  fre¬ 
quency  bandwidth  for  which  (5.10)  holds.  To  examine  this,  ignore  the 
resistive  term  in  (4.2),  and  use  the  normalized  admittances  Bo»  B^  defined 
in  [ 5 ] .  Then 


(5.11) 


where  T  is  the  transmission  coefficient  for  infinite  rib  impedance.  Now 
generally  T  is  small,  while  B  is  then  close  to  its  value  (%)  in  the 

OO  o 

absence  of  fluid  loading,  and  then  it  follows  that  fluid  loading  effects 
on  both  the  magnitude  and  bandwidth  of  the  "resonant"  transmission  are 
negligible.  The  same  applies  at  frequencies  around  coincidence,  where 
T  is  close  to  unity;  in  that  case  we  have  T  ~  1  regardless  of  whether 
the  resonance  condition  ui  =  is  met  in  this  frequency  range  or  not.  If 

we  take  the  very  low  frequency  limit  u  4  0  we  find  that  B  vanishes  like 

U  % 

to  ,  while  we  assume  that  k  varies  like  w  ,  as  it  would  for  a  thin  plate. 

m 

Then  the  second  term  dominates  in  both  numerator  and  denominator  of 
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(5.11),  and  so 


T  ~  T 

OO 

as  to  **■  0.  Thus  for  low  frequencies,  fluid  loading  effects  (represented 
by  the  variation  of  with  to)  again  substantially  "mend"  the  impedance 
discontinuity,  and  except  for  resonant  rib  conditions,  the  low-frequency 
transmission  is  essentially  identical  with  that  for  normal  incidence  and 
infinite  rib  impedance.  At  the  resonant  rib  condition  it  appears  that 
T  ~  1  because  of  mechanical  coupling,  and  fluid  loading  effects  are  not 
able  to  significantly  offset  that  coupling. 

6.  CONCLUSIONS 

It  is  hoped  that  this  analytical  study  of  a  very  simple  model  con¬ 
figuration  will  contribute  to  understanding  of  fluid  loading  effects, 
both  by  providing  insight  into  the  physical  processes,  and  also  by 
complementing  numerical  studies  (such  as  reported  in  [1]  and  [2])  which 
tend  to  run  into  difficulties — even  in  simple  configurations — precisely 
where  fluid  loading  effects  are  most  pronounced  and  most  interesting. 
Highly  accurate  predictions  for  specific  cases  have  not  been  given, 
though  the  analytical  results  given  here  have-  been  put  in  a  form  where 
all  the  sensitive  behavior  has  been  explicitly  extracted,  leaving  simple 
functions  which  can  indeed  be  readily  computed  in  any  specific  case. 
Purely  numerical  computations  of  Fourier  integrals  for  fluid-loaded 
plate  dynamics  can  (indeed,  are  bound  to)  find  great  difficulty  in 
resolving  the  balance  between  structural  and  fluid  coupling;  as  remarked 
above,  this  happened  in  the  work  reported  in  [1]  and  [2],  both  at  low 


frequencies  and  close  to  the  coincidence  condition.  The  perturbation 
method  approach  shows  why  such  difficulties  are  likely  to  arise,  shows 
how  a  balance  is  achieved  between  various  competing  physical  mechanisms 
in  different  frequency  ranges,  and  provides  (at  the  least)  scalings  and 
estimates  of  effects  which  should  be  useful  in  computational  studies  of 
more  realistic  configurations. 

In  Part  II  of  this  work  we  intend  to  pursue  the  application  of  these 
methods  to  structures  with  two  or  more  inhomogeneities. 
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